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Figure 3. Vomiting time versus dose rate. 
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Source: Baranov and Guskova [1988]. 


Figure 4, Vomiting time versus dose. 
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Onset of prodromal symptoms related to dose. 


Accident Brucer 1959; Cronkite et al. 1956: Fanger and Lushbaugh 1967: HUbner and Fry 
1989: Karas and Stanbury 1965: Laumets 1965; Lushbaugh 1973: Thoma and Wald 
1959; Wald and Thoma 1961. 

Therapy Court Brown 1953: Miller et al. 1958: Rider and Hasselback 1968: Salazar 
et al. 1978. 

Composite Gerstner 1958; Cerstner 1970; Glasstone and Dolan 1977: Laumets 1965: Lushbaugh 


Expert Opinion 


1973, NATO 1973; Warren and Grahn 1973; Zellmer 1961. 


Anno 1983, Fanger and Lushbaugh 1967: Gerstner 1958; Glasstone and Dolan 1977: 


Saenger 1963: Thoma and Wald 1959: Wald and Thoma 1961. 


Vomiting onset time (h), referenced from the end of ex- 


posure, was found to be related to dose (cGy) given by, 


where 


T = 10% » pf 
Data Category a A 
Accidents (23 pts.) 2.1830 -0.77260 
Therapy (25 pts.) 1.4611 -0.56271 
Combined (52 pts.) 1.8641 -0.69022 
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Figure 22. Blood count levels following single acute 
radiation exposures. 
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Figure 27. LDEg versus dose rate--three parameter model for mouse. 
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Figure 28. Dey versus dose rate--two parameter model for mouse. 
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Figure 29. LD versus dose rate--two parameter model for mouse, 
from Thompson and Tourtellote (1953) only. 


a oan ae es 0) Ce ee eT | Lot teed nt 


§00- 
LOg, = 190 (1 + (1.27/r)9-94) 
400 
300 
- 
2 
4s 
qQ 
bs] 
200 
% Edmonson snd Batchelor 1971 
a Taviog etal, 1969 
100 - & Pane es at 196K 


Traum (995 


0 TOOT Pep ep tat Lh) Se See en eT ee ee 8 | 
0.1 0 100 1000 10! 
Dose rate (rad’h) 

Source Arrwworth and Leong | 1966) Ainsworth et 2) 11968!. Hankset ai (1966! 
Nachtwey, Aintworih, and Leong 11967], Page (1968), Page otal 11985! 


Figure 30. Dey versus dose rate--three parameter model for sheep. 
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Figure 31. LDco versus dose rate--two parameter mode] for sheep. 


An animal radiation study by Soviet researchers Grigorev, Gorlov, 
and Shafirkin {1978] was translated into English and reviewed. In 
that study, the effects on LD5q of chronic (constant) dose rate ex- 
posure, ranging from 1 to 932 rads/day, was reported for a variety of 
animal species, including mice, rats, guinea pigs, rabbits, donkeys, 
sheep, goats, dogs, and monkeys. Empirical modeling was performed 
that included an extrapolation to man based on a polynomial relation- 
ship linking basal and water metabolism, animal weight, period of 
semirecovery from radiation damage, period of maximum white blood cell 
(WBC) depression, and lifetime of ery*hrocytes. Biological effective- 
ness was defined in terms of a dose-rate dependent coefficient that is 
the ratio of the effective dose D, (i.e., acute LD5g) to the accumu- 


lated dose, D, given as, 


Table 19. Parameters for LD5g versus dose rate. 


Figure 
Species s R Do dgf (RMS/dgf)1/2 No. 
Swine 0.9394 + 0.0321 21.74 229 14 32.51 9 
Swine 1.0 18.77 235 15 32.32 10 
Goat 1.0206 + 0.2160 2.089 240 4 58.94 9 
Goat 1.0 2.146 239 5 52.76 10 
Dog 0.7805 + 0.0913 5.482 263 14 31.78 9 
Dog 1.0 4.074 270 15 33.12 10 
Sheep 0.6887 + 0.1931 1.960 162 8 23.33 9 
Sheep 1.0 1.614 170 9 24.29 10 
Sheep? —s- 0.9423 + 0.1472 1.272 190 4 1721 14 
Sheep@ 1.0 1.258 192 5 10.21 15 
Mice 0.5995 + 0.0856 9.307 873 24 202.1 11 
Mice 1.0 8.717 959 25 202.4 12 
Mice 1.0 6.173 772 4 236.0 13 


4Sheep--without data from Hanks et al. [1966]. 
bMice--data from Thompson and Tourtellotte [1953] only. 
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where r is the dose rate in rads per hour; b = 0.29 for small animals 
(mice, rats, and guinea pigs); and b = 0.22 for large animals and 
humans. When the dose rate, r, is less than or equal to the ap- 
proximate values Listed below, D/Dg can be expressed by the simple 


relationship, 


D/D, =lfar , 


where: 


a(h/rad) 


Mice 0.055 
Rats 0.084 
Large animals 0.168 
Monkeys 0.36 
Man 0.48 


Although our review of the dose-rate dependent relationships 
above does not include comparisons with other empirical models, such 
as the form for LDsg used to fit the animal data shown in Figs. 9 
through 15, and those discussed below (Constant Dose Rate Models). 
Grigorev, Gorlov, and Shafirkin [1978] point out that D/Dp does not 
depend on the dose rate beyond 720 rads/h, as can be seen in the first 


equation above. 


2.4.2 Gastrointestinal Damage. 
Using the split dose technique, Krebs and Leong [1970] performed 


a study with mice to determine the effect of constant exposure rates 


on the gastrointestinal LDs5g9/5 for both 60¢o and 250 kVp X-ray 


irradiation. A plot of their results are given in Fig. 32 that indi- 
cates a factor of about two increase in LDsQ/5 when the dose rate 
decreases from about 8400 to 240 R/h for the 250 kVp X-rays and from 
about 6700 to 52 R/h for 60¢9 gamma rays. The curves, based on probit 
regression fit of mortality, are separated by a factor of 1.48 which 
represents the effective RBE between the two types of photon radiation 
over the range of exposure rate. It should be mentioned that the two 
lowest dose rate points for 60co irradiation 52 and 93 R/h were 
developed from actual exposure periods of only 18 h plus a required 
“topping dose" given at a high dose rate (8400 R/h) to extend the 
radiation to lethality. Based on this data, however, the correspond- 
ing inferred exposure periods would be about 35 and 67 h. It is 
possible that the LDs50/5 values would have been significantly higher 
at those two low dose rates had the exposure periods been allowed to 
full-term lethality. If so, the exposure periods would have been in 


excess of those estimated from their actual 18-h exposures. 
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Figure 32. Effect of dose rate on gastrointestinal LD. 9/5 in mice. 


All actual exposure periods in the Krebs and Leong work were well 
short enough to avoid proliferation of the intestinal epithelium to 
ensure that the LD5o/5 increase with decreasing dose rates only in- 
volved Elkind-type repair. Accordingly, an injury repair model was 
constructed on the basis of the measured recovery kinetics that as- 
sumed exponential repair. With the model, they were able to account 
for the LD59/5 behavior with increasing exposure time and estimate a 
recovery half-time of 23.4 min. 

Krebs and Leong also estimated that some 40 to 50 percent of the 
injury is irreversible. That accounts for an apparent maximum finite 
limit of slightly more than twice the minimum value approached by the 
LDso/5 as the exposure rate becomes low in the absence of other repair 
mechanisms. However, tnis limit may be only transitory due to the 
onset of proliferation which might have been revealed for longer 
exposure periods. For example, as pointed out above, the lowest dose 
rate point (52 r/h) corresponds to an exposure period of 67 h (at 
least) or 2.8 days. In view that the fractionation studies of Withers 


and Elkind [1969] in mice provide clear evidence of proliferation 
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within that period, one would expect a correspondingly significant 
increase in the LDs0/5 associated with that dose rate (i.e., 52 r/h). 
Moreover, Withers [1972] gives evidence of crypt cell proliferation 
for 50 and 60 cGy/h based on constant dose rate exposures. 

Travis et al. [1985] exposed mice to TBI from 60co radiation to 
assess the effect on lethality of constant dose rate, ranging from 60 
to 1500 cGy. Death was scored at ten days to determine the dose rate 
dependence of deaths from the gastrointestinal syndrome. The LDs5q/5 
increased from about 1200 to 2050 cGy (about a factor of 1.7) for a 
decrease in dose rate from the high to the low ends of the range (1500 
to 60 cGy/h). However, the LDsg 5 values are significantly lower by a 
factor of about 1.5 than those given by Krebs and Leong {1970] over 
about the same dose range. The difference could partially be at- 
tributed to different strains of mice used, although deaths scored at 
ten days by Travis et al. rather than within five days compared to 
Krebs and Leong, may have also been affected by damage to the 
hematopoietic system; pathologic findings by Travis et al. indicated 
changes in the jejunum, ileum, and rectum were minimal, and when they 
occurred subjectively, they appeared as a loss of crypts in both the 
jejunum and ileum. Also, in addition to source-to-target distance, 
lead shielding was used to reduce the source radiation to the desired 
exposure rates. This can produce a lower photon energy spectrum from 
scattering, resulting in an increase in RBE that would contribute to 
lower LD5g/5 values. 

Dutreix et al. [1979] described studies of the effect on mice of 
constant dose rate exposures, ranging from 120 to 6000 cGy; they 
indicate a factor of 1.83 increase in LD59/5, somewhat similar to that 
obtained by Krebs and Leong [1970] and Travis et al. [1985]. However, 
the actual LDsg/5 values are again significantly lower than those of 
Krebs and Leong over a corresponding dose range--a factor of about 1.6 
lower at the high dose rate end and 1.8 lower at the low dose rate 
end. Again, some of the difference may be attributable to a different 
strain of mice, particularly in view that the same strain (Balb/c 
mice) used by both Dutreix et al. and Travis et al. were in agreement 


in LD5g values over corresponding dose rates and appeared to be well 


within experimental error. Although not specifically indicated, 


presumably the LDs5g values given by Dutriex et al. involved 60¢o 


irradiation. Whether or not shielding material was employed to at- 
tenuate source radiation exposures, which would change the photon 
spectrum, was not indicated. 

Dutreix et al. [1979] also provided data to demonstrate the 
reduced effectiveness of protracted radiation in terms of a single 
acute dose equivalent of 10 Gy. The data are plotted in Fig. 33 for 
both constant dose rates (connected by the dashed line) and frac- 
tionated exposures (connected by the solid line). The time axis 
includes a three-hour interval which separates dose fractions all 
given in 10 min or less (the three-hour intervals allow adequate time 
for intracellular repair to take place) It is interesting to note 
that over the same period, the continuous dose rate exposure provides 
more protection than fractioned exposures since the acute single dose 


equivalent is a measure of residual damage. Also, the single acute 
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Figure 33. Acute single LD, ys dose equivalent to a 10 Gy protracted 
dose for constant joes rate and fractioned exposures. 


dose equivalent appears to tend to an asymptotic limit with time 
(increasing number of functions) that would appear to be significantly 
higher than that for constant continuous radiation exposure. However, 
it would be difficult to assure this based on the limited data in 

Fig. 33. These data also point out that it would be inappropriate to 
infer the sparing effect of protracted radiation exposure based on 
dose-rate-averaging fractionated exposure. For example, 315 cGy/h 
would correspond to the dase rate obtained by averaging over the 
treatment time for the 2 x 5 Gy fractioned exposure where the acute 
single dose equivalent is larger than that for continuous exposure 
(i.e., 8 Gy compared to 7 Gy). Similarly, the dose rate averaged over 
the treatment time for the 4 x 2.5 Gy fractioned exposure is 109 cGy/h 
where the acute single dose equivalent is 7.5 Gy compared to only 6 Gy 
for continuous exposure. 

Using 200 kVp X-rays, Withers and Elkind [1969] demonstrated 
intestinal radiosensitivity based on their milestone fractionated 
abdominal radiation studies in mice. They developed parallel data for 
the endpoints of lethality and jejunal crypt cell survival and deter- 
mined that 50 percent animal survival within five days corresponds to 
a crypt cell survival fraction of around 0.002. Assuming a random 
crypt cell survival among the crypts, about 1/4 of the crypts will 
survive with proliferating cells. In normal mice, LDs5g/5 (50 percent 
lethality at five days) generally corresponds to about 1/3 of the 
crypts surviving where stem cell survival is 0.003 [Withers, 1989]. 
Potten and Hendry [1983] indicate a higher value of 0.4 (range 0.3 to 
0.5) for fraction of surviving crypts at the LD5g,/5 for mice. 

Fractioned radiation was given at a high dose rate (78 Gy/h) over 
two brief periods (5 to 11 minutes) separated by time intervals rang- 
ing from 2 h to 21 days. The first fraction was 660 rads followed by 
a second one of 1415 rads (total of 2075 rads) for in vivo cell 
colony counting. For lethality assessment, the second fraction was 
given to mice in graded amounts in order to calculate the LDsg 5 from 
probit analyses of data that straddled the LD59/5. The single dose 
LDs09/5 was estimated to 1061 +24 rads, and for split doses, it was 


always higher, as shown in Fig. 34. 
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Source: Withers and Elkind (1969} 
Figure 34, Recovery comparison of LDeg 5 (open circles and 


dashed line) and cell survival (x-symbols at 4-, 
8-, and 24-hours) endpoints; the origin of both 
scales refer to 1060 rads acute single exposure. 


The data points in Fig. 34 (open squares with standard error 
bars) along the dashed "eye-fit" line are the LD5g/5 data that cor- 
respond to the dose fractionation increment, Dj-D,. In terms of 
lethality, it represents the increment of total dose required with 
increasing fractionation interval in order to maintain 50 percent 
lethality in mice at five days. The zero-dose increment value at the 
origin corresponds to a single acute dose of 1060 rads. The X-symbols 
lying dose to the LDs9/5 recovery curve at the 4~, 8~, and 24-h frac- 
tionation intervals are corresponding jejunal crypt cell survival 
ratios (right-hand ordinate) calculated at 1060 rads from expanded 
second-dose fraction cell survival data. The cell survival ratio of 
unity is at zero-dose increment that is normalized to the LD59/5 value 
of 1060 rads for a single acute exposure (cell survival fraction is 
about 0.002). In terms of cell survival, the dose fractionation 


intervals represent the increments of total dose required with in- 


creasing fractionation interval in order to maintain a survival ratio 


of unity (or the cell survival fraction at about 0.002). The 
similarity of the recovery data for cell survival and the LDso ‘5 in 
Fig. 34 is evidence for a causal relationship between 5-day radiation 
death and lethal injury to stem cells of the intestinal mucosa. 
Studies have been done using relatively high dose rates to 
characterize changes associated with intestinal radiation death or to 
determine possible causes. Matsuzawa and Wilson [1965] exposed mice 
to 3000 R of 250 kVp X-rays at a rate of about 40 R/min (-23 Gy/h) and 
determined a mean survival time of 3.5 days ranging from 2.9 to 
3.8 days for conventional mice. For germfree mice, the mean survival 
was 7.3 days with a range of 6.4 to 7.7 days. Based on histological 
examination, they also obtained epithelial cell counts for the crypt 
and villi expressed as percentage versus post-irradiation time shown 
in Fig. 35. The curves show the progressive mucosal denudation with 
time that correlate well with lethality for both germfree and conven- 


tional mice, Based on the progress of thymidine -3H labeled cells, 
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Source: Matsuzawa and Wilson [1965]. 
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Figure 35. Cell counts of mouse ileum after a sterilizing dose of 3000 R. 
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cell transit times were estimated for conventional and germfree mice, 


given below. 


Transit Time (days) 


Conventional Mice Germfree Mice 
Crypt to villi junction 0.5 1.6 
Villi junction to tip 0.75 3.55 
Total (crypt to villi tip) 2.1 4.3 


Cell movement from the base to the tip of the villi was found to 
be linear with time. Also, both in the conventional and germfree 
mice, the normal transit time (i.e., approximate lifetime of the 
intestinal mucosal cells) is about 60 percent of the mean survival 
time. According to Matsuzawa and Wilson, the difference between the 
life span of the cells and time required for denudation (approximately 
mean time to lethality) indicates a change in the life span of the 
cells effected by a decrease in cell population. 

Jackson and Geraci [1986] performed irradiation studies using 
conventional, pseudomonas-contaminated and GI-decontaminated rats to 
investigate the pathophysiological causes of radiation-induced 
gastrointestinal death. They employed fission spectrum neutrons from 
a TRIGA reactor and approximately 8 MeV average energy neutrons from 
an accelerator. Gamma irradiation was also done with !3’cs radiation 
and 60¢o9 radiation at dose rates of about 25 to 26 Gy/h. They-.con- 
cluded that the inability of the denuded mucosa to absorb fluid and 
electrolytes and consequent hypovolemic shock was the major mechanism 
involved in producing intestinal radiation death. 

Figure 36 and Table 20 show some results for LDs9/5 and median 
survival time, respectively, given by Jackson and Geraci. The RBE of 
cyclotron neutrons appears to be about two for gastrointestinal 
lethality regardless of treatment condition; for fission neutrons, the 
RBE is almost three (2.8). However, no significant difference was 
found in the median survival time for conventional rats regardless of 
the radiation source including gamma rays from either 137cs or 690 
and neutrons from either the cyclotron or TRIGA. Also, the time of 
death is not explicitly dose dependent but given a sufficient amount 


of gastrointestinal damage (i.e., lethality threshold), the time of 
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Figure 36, 


Survival time of untreated, GI-decontaminated, and 
Pseudomonas~- infected animals after exposure to various 
doses of Cs-137 gamma rays and cyclotron neutrons (each 
point represents the median survival time of 8 to 40 rats). 
Source: Jackson and Geraci, 1986. 
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Table 20. LDs5o/5 day of pseudomonas-infected GI, conventional, 
and GI~-decontaminated rats. 
Source: Jackson and Geraci, 1986 


LDsosday 
Treatment Radiation (Sy) 

Pseudomonas 137Cs y rays 13.7 (13.1-14.3)4 
Pseudomonas C-neutrons 5.85 (5.11-6.69) 
Conventional 137Cs y rays 14.1 (13.4-14.7) 
Conventional C-neutrons 6.25 (6.10-6.40) 
Conventional F-neutrons© 5.0 (4.72-5.30) 
Decontaminated 137Cs y rays 15.0 (13.9-16.2) 
Decontaminated C-neutrons 7.40 (6.79-8.09) 


avValues in parentheses are 95 percent confidence limits. 
C-neutrons = Cyclotron produced neutrons. 
©F-neutrons = TRIGA produced fission neutrons (1.2 MeV mean energy). 


death proceeds according to degenerating physiological processes. 
When radiation exposure is protracted, the damage threshold is offset 
by recovery mechanisms which may either prevent or delay lethality. 

Very little real data exist for lethality in humans that can be 
attributed purely to gastrointestinal injury. In fact, aside from 
being a predominate cause, it probably can never be expected to be the 
sole cause of lethality following TBI owing to other accompanying body 
injury. The most comprehensive set of lethality data for humans 
involving gastrointestinal injury is from the Chernobyl accident 
victims listed in Table 21 given by Guskova et al. [1988]. Dose 
estimates range from 10 to 12.5 Gy and the time of death from 10 to 
20 days following the accident; the dose rates are indicated to be at 
a high level of at least 2 Gy/h or more {Baranov and Guskova, 1988]. 
Rather severe injury to the skin from beta radiation was a competing 
lethal effect as indicated in Table 21. The fact that some 
lethalities extended up to 20 days might reflect the extensive medical 
care given. 

Table 22 gives estimates of dose and time of gastrointestinal- 
syndrome-lethality for man based on some pre-Chernobyl experience 


(with the exception of that indicated by Guskova et al., 1987). The 


Table 21. Mortality of Chernobyl accident vi: cims 
with skin and intestinal injuries. 


Case Dose (Gy) Day of 
Number (Marrow) Treatment? Death 
3 12.0 BMT 17 
4 11.8 BMT 18 
10 11.1 FLT 14 
14 10.9 FLT 18 
15 >10.0 FLT 14 
17> 10.0 BMT 18 
20 12.4 FLT 17 
23 13.7 FLT 15 
26> 12.5 20 

2097 (Kiev) 10.2 10 


4BMT = bone marrow transplantation 
FLT = fetal liver (cell) transplantation 
involved mycobacterial sepsis 


b 
accident cases given by Fanger and Lushbaugh [1967] are thought to be 

predominately due to gastrointestinal injury. Even though some of the 
dose ranges in Table 22 extend to 7000 cGy, Fanger and Lushbaugh 


[1967] give evidence for cardiovascular shock as the radiation 


Table 22. Gastrointestinal syndrome lethality in man. 


Dose (cGy) Postexposure Time (days) Source 
1114, 1910 9, 10 Fanger and Lushbaugh (1967) 
1000-2000 8-165 Guskova et al. (1987) 

21000 <14 Maisin et al. (1971) © 

>1000 7-10 Fajardo (1982) 

670-67004 7 Ingram (1969) 

750-2000 7-14 Anno et al. (1989) 

670-67004 6-10 Lushbaugh (1973) 


“Two accident cases cited from Hemplemann et al. (1952) and 
Kurshackov (1962). 

bIntestinal changes most apparent, 8-12 days; total loss of 
epithelium, 10-16 days. 

©Based on histologic examinations from autopsies performed on eight 
Hiroshima and Nagasaki bombing victims (substantial dose uncertainty) 

dBased on a 0.67 internal-body-dose conversion factor for 
1000-10,000 rads. 
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syndrome causing death based on pathologic changes in two accident 


victims who received doses of 4500 and 8800 cGy; furthermore, the 
times of death, 35 and 49 h, respectively, following exposure, are 
much too early for pure gastrointestinal injury to manifest in 
lethality. 

We conclude that lethality solely from gastrointestinal injury in 
man exposed to TBI may exist in a limited range, perhaps ~8 to 12.5 Gy 
equivalent acute dose. Lethalities can be expected to occur within 
one to two weeks. However, GI injury may play an important part in 
lethality for doses on either side of that range. Also, when radia- 
tion exposure is protracted, the recovery potential of the intestinal 
mucosa can produce an increase in the total dose level required to 


effect lethality. 


2.5 RADIATION INJURY AND RECOVERY. 

Radiation injury and recovery research in mammalian species 
relevant to this effort is primarily based on the lethality endpoint 
from radiation-induced bone marrow cytopenia. The empirical tech- 
niques for the measurements required have been rather well established 
in a number of mammalian species. In other tissues, such as the 
intestines, in vivo injury and recovery assessment techniques with 
regard to the effect manifested in the whole organism have not been as 
well established. Accordingly, there exists a much larger body of 
literature with regard to demonstrating injury and recovery due to 
hematopoietic effects. However, compared to the hematopoietic system, 
the much earlier and more rapid rate of recovery of the intestines 
from radiation injury has been well established based on the work of 
many researchers including Quastler [1959], Quastler et al. [1959], 
Hornsey and Vatistas [1963], Lesher [1967], Withers and Elkind [1968, 
1969, 1970], and Potten and Hendry [1975]. A more detailed discussion 
of intestinal injury and recovery accompanies the gut injury model we 
developed which is described in Section 5. 

The postirradiation injury state, i.e., residual injury or injury 
that remains unrepaired, of an organism depends on the following 


considerations: 
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1. The extent to which the injury is dose dependent 
(i.e., radiosensitivity). 

2. The physiological repair and recovery which can 
counteract the effects of injury. 

3. Recovery must occur within a definite time or 
to a certain extent to effect survival. 

4. Radiation exposure over a period of time induces 
fewer injuries per unit of time, permits repair 
processes to begin earlier, and enables prolifera- 


tion to take place more rapidly. 


Based upon lethality in mice from damage to the blood system, Blair 
[1952] developed a theory that radiation injury is repaired at a 
constant percentage per unit time of the net recoverable injury, and a 
small fraction of this injury is irreparable. This can be expressed 


in terms of net injury (acute single dose equivalent), De as 


2 7 “Bt 
De = DS [f+ (1 - fje "7] 


where f = irreparable fraction of injury, 
8 = recovery rate in percent per day, 
t = number of days, 


Do = single acute dose. 


Davidson [1957] applied this equation to the results of a split- 
dose recovery study in mice by Patterson, Gilbert, and Mathews 
[1947]. He then derived the graph depicted in Fig. 37 where net or 
residual injury is expressed in roentgens. Davidson estimated 
10 percent irreparable injury, and determined that 2.5 percent 
recovery per day would yield a curve that would best fit the data of 
Patterson, Gilbert, and Mathews. He then proceeded to estimate WBC 
recovery half-life for several mammalian species and correlated these 
data with minimum white cell counts in response to an LDspg exposure. 


He found an approximate linear correlation for the mouse, rat, dog, 
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@ = experimental results 
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f 
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10 percent irreparable injury 
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Source: Davidson [1957] . 


Figure 37, Theoretical injury recovery curve with fraction of 
injury assumed unrecoverable. 


and burro, Using human data available to him for WBC counts [Oughter- 
son and Warren, 1956], Davidson estimated a recovery half-life for man 
of 690 h. 

Subsequently, several investigators conducted a number of split- 
dose recovery and hematological studies to determine the validity of 
postirradiation exponential recovery and irreparable injury (Holloway 
et al., 1968; Ainsworth and Leong, 1966; Page et al., 1971; and 
Nachtway, Ainsworth, and Leong, 1967]. Irreparable or residual injury 
was confirmed; it appears to be dose, dose rate, and species dependent 
{Baum and Alpen, 1959; Baum, 1967]. Two examples of residual injury 
in the erythropoietic system are shown in Figs. 38 and 39. Rats were 
repeatedly irradiated with 300 or 400 R X-rays--five times at 90-day 
intervals. After each exposure, red cell production decreased during 
the recovery period. This decrease was significantly greater in rats 
subjected to 400 R X-radiation (Fig. 38). Figure 39 indicates that 
enthropoiesis significantly increases over the first five days follow- 


ing exposure. 
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In Figs. 40, 41, and 42, data are plotted of split-dose recovery 


Figure 38. Maximum Fe? incorporation for 
male and female rats exnosed 
periodically to 300 r of X-rays 
and for female rats irradiated 
with 400 r (isotcne was injected 
5 days postirradiation, 3-month 
interval between exposures). 
Source: Baum and Alpen, 1959. 


studies in hamsters and sheep performed at the Naval Radiological 
Defense Laboratory (NRDL) and rhesus monkeys at the Air Force School 
of Aerospace Medicine Laboratory [Holloway et al., 1968; Ainsworth and 
1967]. 
tial postirradiation recovery, over varying numbers of days, which is 


interrupted by a period of greater radiosensitivity (increased injury) 


ae et a 
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Figure 39. Residual injury of 


erythropoietic system 
in irradiated rats. 
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Figure 40. Hamster--acute exposure. 
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Figure 41. Rhesus monkey--acute exposure. 
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Figure 42, Sheep--acute exposure, 


before final return to recovery occurs. None of these three species 
of animals appears to return completely to preirradiation levels. 
Similar reversals in recovery have been reported for the mouse and the 
rabbit (Leong, Wisecup, and Grisham, 1964]. 

Experiments designed to determine the pathophysiological 
responses of the blood-cell-forming system show postirradiation oscil- 
latory recovery, which might be at least partly responsible for the 
return to hypersensitivity in the above-described split-dose recovery 
studies [Baum, 1967; Baum and Wyant, 1970; and Morley, King-Smith, and 
Stohlman, 1970]. Figures 43 (for the rat) and 44 (for the dog) 
demonstrate these oscillations for postirradiation red cell recovery. 
The finding that the above animals become radiosensitive after a 
period of substantial recovery suggests that a critical organ system 
(perhaps the bone marrow) has undergone some alterations such that 
little additional radiation is required to kill the animal. This, of 


course, questions the validity of Blair's theory of an exponential 
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Source: Baum {1967}. 


Fiqure 42, Residual injury of erythrocyte stem cells 
in irradiated rats--1l7-day observation 


period. 
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Figure 44, Erythropoiesis in normal and in polycythemic dogs 
exposed to 150 rads mixed gamma-neutron radiation. 
Source: Baum and Wyant, 1970. 
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recovery. Leong, Wisecup, and Grisham {1964] report a return to 
radiosensitivity for rabbits about three weeks after irradiation. 
Prior to that, continuous recovery of almost 60 percent was recorded. 
Since at the time of the reversed sensitivity, peripheral leukocytes 
were in the recovery phase and continued to recover, the authors 
questioned an involvement of the hematopoietic system with the return 
to radiation hypersensitivity. However, the postirradiation con- 
centration of leukocytes in the circulatory blood system merely indi- 
cates that the perturbed bone marrow shunts newly formed leukocytes 
rapidly to the circulation for transportation to the tissues to 
protect against invading antigens. It does not necessarily mean that 
normal storage areas have been replenished and that sufficient numbers 
of white cells are available to counteract the deleterious effects of 
an additional dose of radiation. 

In two other experiments (Figs. 45 and 46) conducted at NRDL 
using the dog and the swine, return to an increased injury condition 
is apparently not indicated. However, it must be emphasized that both 
studies were only carried out to day 20 postirradiation, and reversal 
of injury is seen in the sheep, another large mammal, between days 20 
and 30. 

For the swine, an experimental point is available for day 60 
which could not be statistically evaluated. However, the swine may 
have a remarkable capability for recovery and be qu te radio- 
resistant. Figure 47 shows that swine exposed to an acute dose of 
155 eGy given in about 0.5 h have sustained injury at the conclusion 
of the exposure equivalent to 155 cGy [Ainsworth et al., 1968]. If 
the same dose is given in 61.5 h, the remaining injury is only 
64 cGy. Doubling or tripling the dose increases the injury only 
slightly. One possible interpretation of this Finding is that swine 
develop increased resistance to radiation while being irradiated. 
Another may well be that animals sustain the greater part of their 
injury during the early radiation period. Experiments with sheep by 
Still et al., (1969) and Jones and Krebs [1970] seem to give credence 


to the latter interpretation. 
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Figure 45. Qog--acute exposure. 
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Figure 46. Swine--acute exposure. 
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TIME (HOURS) 


Figure 47, Swine--injury and recovery versus exposure time 
at fixed dose rate. 


As indicated in Table 23, Still et al. [1969] gave animals 95 cGy 
at an acute dose rate of 311 cGy/h. Immediately after the completion 
of this conditioning (initial) dose, irradiation was continued at a 
reduced dose rate of 2.2 cGy/h until an LDs59 of 199 cGy was obtained; 
this value was not significantly different from the LD5o/69 of 192 cGy 
observed when sheep received an acute dose rate at 275 cGy/h. 

However, 192 cGy did differ significantly from the LDs5q/gq of 302 cGy 
obtained with a protracted dose rate of 2.2 cGy/h (Ainsworth et al., 

1968]. The hottom half of Table 6 shows similar results observed in 

studies by Jones and Krebs {1970}. 

Figure 48 shows postirradiation recovery of five large mammalian 
species. The dog, sheep, and swine show similar recovery patterns 
over the first 20 days postirradiation, but otherwise, these curves 
are quite dissimilar. It is difficult to extrapolate a recovery 
pattern for humans from these curves unless other response parameters 


measurable in humans exposed to ionizing radiation can be correlated 
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to them; reevaluation of postirradiation leukocyte responses in man 
might be helpful. 

Sheep are the only large animals for which data are available to 
compare the radiation effects of various low dose rates with an acute 
one. Such studies were conducted at NRDL, and the results are 
depicted in Figs. 49 to 53. As may be seen from Table 24, the first 
five conditioning doses employed were similar. It appears that when 
sheep are exposed at below a dose rate of 1 cGy/h, no apparent injury 
remains at the end of irradiation. Above a dose rate of 1 cGy/h, 
injury can be determined at the end of the radiation exposure, which 
appears to increase with increasing dose rate. As the conditioning 
dose increases (186 cGy in Table 24), the remaining injury increases 
in direct proportion. When the conditioning dose is 101 cGy, injury 
accumulates at 0.63 cGy/cGy. However, when the total dose is in- 
creased to 186 cGy at the same dose rate, injury accumulates at only 
0.5 cGy/cGy. Maximum over-recovery, given a similar size of con- 
ditioning dose, increases with dose rate and appears to be propor- 
tional to remaining injury, at least over the dose rate ranges indi- 
cated in Table 24, It appears that the severity of initial injury may 
be related to a humoral or hormonal release that stimulates increased 
cellular production in recovering bone marrow. 

It has been established that recovery and repair are important 
factors that enable mammals to sustain increased radiation injuries 
with protracted exposures. Split dose experiments permit the measure- 
ment of recovery based on the increase in the radiation dose to 
satisfy a specific endpoint (e.g., 50 percent lethality). More com- 
plete measurements obtained from split-dose experiments could 
facilitate the development of relationships for dose protraction in 
animals and eventually extrapolations to humans. If physiological 
processes cause a deviation from a simple recovery process such as a 
return to hyper-radiosensitivity after a period of normal recovery 
(i.e., an undulating pattern), the underlying pathophysiology of the 
system which induces such changes must be considered. Cell cycling 
sensitivity to irradiation exposure, as demonstrated by Bedford et al. 


{1980], may also play a part, which should also be among the underly- 
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Figure 48. Species comparison for acute exposure. 
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Figure 49, Sheep--acute and protracted dose 
(0.63 and 0.32 cGy/h). 
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Figure 50. Sheep--acute and protracted dose 
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Sheep--acute and protracted dose over 
1.8 days (2.4 cGy/h). 
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Figure 52. Sheep--acute and protracted dose 
over 3.2 days (2.4 cGy/h). 
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Figure 53. Sheep--acute and protracted dose comparisons. 
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Table 24. Remaining injury and maximum over-recovery 
in sheep subjected to radiation at different 


rates. 

Conditioning Dose Remaining Maximum 
Radiation Dose Rate Exposure Time Injury Over-recovery 
(cGy) (cGy/h) (h) = (days) (cGy) (cGy) 

108 275 0.4 _ 108 -44 
101 0.32 316 13 0 -19 
101 0.60 168 7.2 0 -26 
OL 1.1 92 3.8 46 -64 
101 2.4 42 1.8 64 -188 
186 2.4 77.5 3.2 93 -75 


ing casual factors considered. 


2.5.1 Constant Dose Rate Models. 

Accumulated injury models were reviewed by Anno and Baum [1986] 
in terms of LDsg endpoint as a function of constant dose rate level 
such as depicted in Fig. 54. Figure 54 also summarizes some of the 
findings of Krebs and Jones {1975] from their review of animal data 
(sheep, dogs, swine, goat, and mice) and suggestions for modeling LDs5q 
response in animals, including extrapolation to humans. They 
developed a relationship LD5q (D59q in Fig. 54) as a function of dose 
rate r, which includes linear and exponential relationships, to model 
protracted radiation response. In their analysis of the animal data, 


Krebs and Jones found, in part, that 


@ lethal dose becomes dependent upon dose rate when the 
time required to deliver it is longer than about 
30 min (otherwise, the LDsg-dose-rate relationship 


is flat where dose rates can be considered “acute"), 
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Source: Krebs and Jones (1975]. 


Figure 54. LD... versus dose rate. 
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@ the LDs5q-dose-rate relationship is linear between about 
20 to 30 rads/h and the acute dose rate level, 

@ between about 2.5 rads/h and 20 to 30 rads/h, there is 
a transition from a linear LDsog-dose-rate relation- 
ship to an exponential one, 

@ less than about 2.5 rads/h, the LD5g depends only on the 
average daily dose rate which can be averaged out over as 


much as two weeks. 


Anno and Baum [1986] compared accumulated human injury models by 
plotting LD59 (rads) versus dose rate r (in rads per h). That com- 
parison, shown in Fig. 55, required some algebraic adjustments to 
obtain the appropriate plotting forms that are given below along with 


the parameters. 
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Figure 55, LDen dose rate models. 


Strandqvist [1944]. 


_ ny [1/0-b)} -[b/(1-b) } 
Deg = Dy + (168r) 


where Dy = 345 rads (acute LDsg) and b = 0,26, 0.52. 


Bateman [(1968}. 


where Dg = 300 rads (fixed), K = 1.64 (rads/h)1/3 based on the 1964 


Mexican accident [Martinez et al., 1964]; and Do = 236 rads, K = 2.29, 
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also based on the 1964 Mexican accident, but with the acute dose 


anchored at Dy = 300 rads for r = 600 rads/h. 


Equivalent Residual Dose [Blair, 1952; Davidson, 1957]. 


Dey = : {exo - oi Ls exp (-ADs9/t)] , 


where ERD (equivalent residual dose) = 300 rads, f (irreparable frac- 
tion) = 0.1, and the repair constant 6 = 0.00104 h7! and 0.002083h7!; 
(for B = 2.5 and 5.0 percent per day, respectively). 


Operational Equivalent Dose, {Home Office Scientific Research and 


Development Branch, 1985]. 


where OED (operational equivalent dose) = D5g (acute) = 450 rads. 

Data from other sources are also individually plotted in 
Fig. 55. The two box-shaped "R42" values are based on the LD59 values 
given in the "Penalty Table" by the National Committee on Radiation 
Protection and Measurements [1974]. For one-week exposure, an LDs5q of 
300 rads (450 R) is given; that corresponds to an average dose rate of 
about 1.77 rads/h (2.68 R/h). For one-month exposure, an LD5g value 
of 400 rads (600 R) is given; that corresponds to an average dose rate 
of about 0.55 rads/h (0.82 R/h). 

The two values marked "BIR" in Fig. 55 are based on information 
from the British Institute of Radiology (BIR) as quoted on p. 84 in 
the British Medical Association report of 1983. 

The rectangles in Fig. 55 reflect dose and dose rate uncertainty 
and are based on the 1964 Mexican accident involving cobalt-60 y-ray 
radiation exposure of five family members [Martinez et al., 1964] that 
resulted in four deaths (+) and one survivor (S). Accordingly, these 


are not LDsg data, but are shown for reference only. 
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Plots of the protracted radiation response models (Fig. 55) show 
a considerable variation in accumulated lethal exposure dose versus 
dose rate. However, with the exception of the Bateman model, those 
plotted suggest a marked increase in LD5g, commencing with dose rates 
less than about 3 to 10 rads/h; the even more rapid increase in LDs5go 
for dose rates from about 1 to 3 rads/h probably reflects cell 
proliferation. Anno and Baum {1986] provide a detailed discussion of 
the model plots. 

A brief review of some suggested lethality-endpoint-based models 
of protracted radiation response illustrates the need for additional 
investigation. Lethality is only one of the endpoint responses of 
interest in casualty considerations: however, the models do predict 
various degrees of biological recovery. Any model selected should be 
better substantiated by more in-depth analysis of available data from 
animal studies and preferably human experience; for example, the 
Goiania, Brazil accident involving protracted exposure to cesium~137 
y-ray radiation from an abandoned teletherapy unit [International 
Atomic Energy Agency, 1988] could yield more clinical information. At 
present, data on arbitrary exposure periods and/or varying dose rates 
are scarce or limited in scope. Consequently, our comparisons of the 
protracted radiation exposure models are based on continuous and 
constant exposure rate levels. However, as Krebs and Jones [1975] 
imply, when average daily dose rates are less than about 2.6 rads/h 
for about 62 rads/day), the exposure history for the 24-h period is 
largely irrelevant. 

Our investigation of the kinds of models reviewed here indicates 
that they cannot be generally applied to other endpoints such as 
prodromal responses without possibly data-supported modifications--and 
then only for constant dose rate exposure. Hovever, such models have 
been applied for prodromal symptomatology of protracted radiation 
involving fallout effects analysis [Knapp, 1965; and Schmidt, 1981]. 
Because the kind of biological recovery illustrated in this review may 
not adequately model other processes (such as a physiological clearing 


action and recovery), a different type of modeling approach is 


generally necessary to accommodate prodromal responses to protracted 


radiation. For example, for dose rates in the therapy range of about 


1 to 30 rads/min (60 to 1800 rads/h), there are indications that 
nausea and vomiting depend more on the total accumulated dose than on 
the dose rate [Baum et al., 1984]. Our assessment of the existing 
models of protracted radiation based on lethality as the endpoint 
reveals that additional study of available data is needed before 
proceeding to a systems analysis approach that is used in military 


operations and planning. 


2.5.2 Operational Equivalent Dose as Injury Accumulation. 

The residual injury accumulation based on the OED formula recom- 
mended for appiication in the UK [Home Office Scientific Research and 
Development Branch, 1985} is illustrated in Fig. 56. The OED formula 
is given by 


OED = D - 150 - 10t (rads) 
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Figure 56. OED (injury, for constant daily 
exposure rate. 
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where D refers to the accumulated dose to the bone marrow, t is the 


exposure time in days, 150 (rads) is a rapid (Elkind) repair or 
recovery value that takes place within the first day following ex- 
posure, and 10 rads/day represents a daily recovery rate. Given a 
dose rate R, and letting D = Rt, the OED (or accumulated injury) can 


be expressed as 


OED = (R - 10)t - 150 (rads) . 


Plots of the OED relationship above for dose rates, R = 10, 20, 
and 30 rads/day are shown in Fig. 56. Since negative values are 
ignored in applying the OED formula, the residual injury would remain 
at zero until it accumulates linearly after specific exposure times, 
depending on the dose rate, as shown. 

The OED was also compared to data derived from experiments with 
sheep {Hanks et al., 1966]. Figure 57 illustrates the comparison for 
radiation exposure protracted over 3.7 days at the rate of 1.1 rads/h 
ox 26.4 rads/day. The initial portion of the bottom curve (large 
dashes, marked OED) is a plot of the OED relationship given above over 
the period of exposure. For postirradiation times where R = 0, it is 
assumed that recovery continues at the rate of 10 rads/day. 

It is clear from Fig. 57 that the OED considerably overestimates 
recovery from injury compared to the data drawn from experiments with 
sheep. The overestimate is primarily due to the 150 rads that purpor- 
tedly accounts for "rapid recovery." When the 150 rads of repair 
recovery is ignored, given by the curve marked OED + 150, the agree- 
ment with the experimental data is vastly improved out to about 
12 days postirradiation. This underscores the most significant 
problem with the OED formula, although the recovery rate of 
10 rads/day appears to be within reason. 

Our review and analysis of the literature provided us with a 
comprehensive assessment of radiobiological injury and recovery 
relevant to acute and protracted radiation exposure. In the section 
that follows, we discuss aspects of the review key to the development 


of modeling approaches for upper and lower gastrointestinal distress. 
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APPENDIX A 
UGIDM EQUATIONS FOR ACUTE AND CONSTANT DOSE RATE EXPOSURES 


This appendix derives solutions of the UGIDM equations for two 


specific cases of radiation exposure: acute (or prompt) and constant 


dose rate. The model variables and parameters are: 


Variables 


depletable reservoir (target tissue) level, Gy 


P = potential toxin, Gy 
A = active toxin, Gy 
R = dose rate, Gy/h71 
D = dose, Gy 
Parameters 
a = potential toxin conversion rate/active toxin production 


rate, hol 
8 = active toxin clearing rate, nol 
uw = depletable reservoir reconstitution rate, h7l 
Cy = initial reservoir level, Gy 
Do = characteristic target tissue dose (=C,), Gy 
Ago,5 = half-maximum value of A, Gy 


y = severity response slope 


A.1 ACUTE EXPOSURE. 


The UGIDM differential equations for acute exposure are, 


depletable reservoir: a6 »(c,-C) 


tential toxin: _ = =e 
potentia oxin ac aP 
active toxin: a = aP-fA . 
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(A.1) 


(A.2) 


(A.3) 


Following an acute exposure, the initial amount of potential toxin 


derives from the depletable reservoir in the UGIDM. 


This can be 


determined by taking the appropriate limits of the solution to the 


differential equation for the depletable reservoir at a constant dose 


rate, R, given by, 

R 
a" >, C+ #(c,-¢] ‘ 
or —— + KC = ue, , Kes R/D, + 


and then, 4 (c_Kt] = pc et 


Ce t 
Integrating, | a(ceX*) = BC, | eft ae , 
° 


and solving for C(t) gives, 


pC 


ae (a -e kK} + ear 


C(t) = K 


Taking limits of the dose rate (R) and time (t) along the fixed 


relationship, D = Rt (where R/Dg>>y), 


fim C(t) = coe Po 


R 7 @ 
t 7-0 


The initial amount of potential toxin Py is then, 
P =C - ceo : 


fe] ° 


The solution to Eq (4.2) is, 
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(A.4) 


(A.5) 


(A.6) 


P(t) = bes : (A.7) 


For the active toxin, Eq.(A.3) is, 


dA 
dt + BA =aP , 
d pt -at Bt 
and, S (re } = aP ie ee : 
where ePt is the integrating factor. 
ga* t 
Integrating, a{ae?*] = oP, | e(P-a)ta. 
A 0 
fe) 


and solving for A(t) gives, 


aP 


A(t) * ie (.-at_.“Ft] (A.8) 


The maximum UG severity level following an acute exposure is obtained 
by maximizing A(t). That is, for dA/dt = 0, the time that A is a 


maximum following acute exposure is obtained as, 


t = [zn 4 / (Bra) . (A.9) 


max 


Then the maximum UG severity is, 


Saooe Ia {1 -exp|-[A(taax) / os] "}} (A.10) 


A.2 CONSTANT DOSE RATE EXPOSURE. 
The UGIDM differential equations for constant dose rate exposure 


are, 


_,. dC _fR _ 
depletable reservoir: Ae = | C+ ph (S, c) (A.11) 
; at dP R),_ 
potential toxin: ac. = >, Cc aP (A.12) 
dé 
active toxin: ae aP - BA (A.13) 


The solution to Eq. (A.11) is given by Eq. (A.5) above. Multiplying 
by the integrating factor e®¢, Eq. (A.12) may be given by, 


Integrating, 


at 
Pe t uC 
| a(pe"] take ae ( = okt} +C oR tae . 
D K o 
P ° 
re) 
and solving for P(t) gives, 
RC 


P(t) = = [E [a - e*) - E) (-**-2°*5}] + Poe. (A.14) 
°o 


Multiplying by the integrating factor eft, Eq (A.13) may be given by, 


a-p 


Sott) = (A fal em) - fx] ote geet} 
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Integrating, 
Bt ‘ 
(ae**] = aP | elFrOE ae 
fe) 


e 


A 
° 


ct 


ree) {,[eF : Fe] : o() [260 (Fa)e]) ae 


and solving for A(t) gives, 


wo HE BO eo ego 
Oo 


Ke ee aP a _ = 
“he EJ] FEO cae as 


Eq. (A.15) gives the active toxin value as a function of exposure 
time. When Py, Ap ~0, they would be residual values from previous 
exposure(s) referenced from t=0 at the beginning of the current con- 
stant dose rate exposure; otherwise, P,, Ag = 0 and the last two terms 
are zero. For this condition, the time that A is a maximum during 
irradiation (which may the time irradiation ceases) can be determined 
for dA(t)/dt = 0; then when Py, A, = 0, a transcendental relationship 
results which must be solved for t = tmayx numerically, i.e., of the 


form below, 
ao (Kr) . Ae (278)t = 8B 


where, the A, B, and K are constants that depend on the dose rate R 
and model parameters. 


The active toxin A(t’) after irradiation ceases is, 


aP’ ' -~Aarr = ’ 
A(t’) = re Fae ccs ke Are a (A.16) 
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where, t’' is the time measured after irradiation ceases, starting from 
zero at that point, and Pg and Ag are the existing values then for P 
and A; they are given by Eqs. (A.14) and (A.15), respectively, for Po, 
A, = 0. The time A(t‘) is a maximum after irradiation ceases is 


determined for, dA(t’)/dt’ = 0, given by, 


A‘ 
pts (hea) mf - 3 és} (A.17) 


oO 


The maximum severity level of UG distress is then given by S(t‘ may) 


according to Eq. (A.10). 
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APPENDIX B 
INCIDENCE OF NAUSEA AND VOMITING WITHIN 24 HOURS 
IN NUCLEAR ACCIDENT EXPOSURES 


A log likelihood analysis was performed based on 40 different 
cases of accidental acute exposure of humans to nuclear radiation. 
Five different models including the normal, lognormal, Weibull, logis- 
tic, and loglogistic forms were applied to fit nausea and vomiting 
response data for the first 24-h postexposure. Two sets of calcula- 
tions were performed based on somewhat different internal whole-body 
dose estimates for 13 of the 40 cases. Most of the accidents involved 
mixed neutron and gamma radiation exposure, and an RBE = 1 was assumed 
for analyzing the nausea and vomiting responses. 

ED5g estimates range from 148 to 177 cGy for nausea and 160 to 
190 cGy for vomiting when all five models and both data sets are 
considered, All of the models fit the data reasonably well; the 
lowest x2 values were obtained for the lognormal, Weibull, and log- 
logistic models. Based on data set I and the lognormal model, EDso 
values of 157 (+2§) cGy were obtained for nausea and 170 (+88) cGy for 
vomiting. The EDsg values for the alternative data set are essen- 
tially the same, differing by about 9 cGy. 

Calculations were performed to estimate the incidence of nausea 
and vomiting for acute radiation exposure in order to establish an 
acute (or high dose rate) basis relevant to our study of protracted 
radiation exposure. The population we chose to analyze for this 
purpose consists of 40 nuclear accident cases documented in the 
literature where total body radiation exposures were all delivered in 
the order of minutes or less at rates exceeding ~2000 cGy h71, Much of 
the large body of clinical information for total body irradiation 
(TBI) of patients with various malignant diseases also involves acute 
radiation exposure rates. However, the prodromal response is altered 
to varying degrees due to chemotherapy often administered prior to 
TBI, premedication (antiemetics, sedatives, analgesics, and steroids) 


as well as the health condition and well-being of the patients. 
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Accordingly, we avoided the use of clinical TBI data because of these 


limitations. 


B.1 DATA. 

Nuclear accident data used for this analysis of the incidence are 
given in Table 31. The data covering a variety of incidents are 
adapted from the literature. Each of the 40 cases are for an in- 
dividual accidentally exposed to acute ionizing radiation and are 
labeled (column 1) designating geographical location and specific 
clinical case (see Table 31 footnote); in some instances, more than 
one individual was exposed during a given incident. Specifically for 
the regression analysis, we utilized the dose values given in column 4 
and the response data in columns 5 and 6 where "X" indicates a posi-~ 
tive quantal response. 

Most of the case data in Table 31 were originally compiled by 
Thoma and Wald [1959]. However, based on subsequent dosimetric 
analysis, some of the originally reported dose values were revised 
downward to reflect internal dose estimates to the midline of the body 
or the mean bone marrow dose. Table 31 includes the revised dose 
values as well as the doses for additional accident cases given in the 
other references [Lushbaugh, 1969; HUbner and Fry, 1980; Klener et 
al., 1986; and Fanger and Lushbaugh, 1967] listed in Table 31. More 
recently, Baverstock and Ash [1983] also performed additional 
dosimetric analysis of the Oak Ridge (OR) and Yugoslavian (Y) 
accidents. They estimated somewhat lower dose values for 13 of the 40 
accident cases, as given by the values in parenthesis in Table 31. 

Our analysis considered both sets of dose data. 

A large majority (about 80%) of the accident cases shown in 
Table 31 involve substantially varying proportions of mixed neutron 
and gamma radiation (column 3) due to nuclear criticality incidents. 
Data are lacking, however, to assign any neutron RBE effect to the 
nausea and vomiting response in humans. Accordingly, we have tacitly 


assumed an RBE of unity for those prodromal responses, and the neutron 


and gamma doses (in gray units) are utilized additively. 


Table 31. 


Nuclear accident cases--nausea and vomiting within 24 hours, 


Dosimetry 
Internal Prodromal Response 
Case* Reference? n/y Ratio Dose (cGy)* Nausea Voaiting 
La2 1,2 a0 8.1 -- -- 
LALO 1,2 3.5 9.0 -- -- 
Aa 1 0.08 10.8 -- 2s 
Lag 1,2 3.0 12.0 -- -- 
Las 1,2 3.0 16.0 -- -- 
OR8(H) 1.2.3 0.36 22.8 se as 
La? 1,2 3.7 42.0 -- ag 
aj 1 o.1 60.5 -- “7 
LA6 1,2 4.6 62.0 -- -- 
OR6(F) 1,2,3.4 0.36 68.5(66.5)4 -- -- 
OR7{6) 1,2,3,4 0.36 68.5(66.5) -- -- 
P(A) 2 0.0 100.0 -- ts 
A2 3 0.083 125.6 -- -- 
c 5 0.0 140.0 -- -- 
Y6(B) 1,2,4,6 0.28 145.0(127.0) x -- 
Al 1 0.098 159.2 x x 
UT/CARL 2 0.0 165.0¢ -- x 
Las 1,2 6.4 192.0 x x 
NJ(2) 2 0.0 200.0 x x 
YS(H) 1,2,4,6 0,26 226 .0(201 .0) Xx X 
ORS(E) 1.2.3.4 0.36 236 .0(225.0) x -- 
OR4(B) 1.2,3.4 0.36 270.0(265.0) x x 
Y4(G) 1,2,4,6 0.28 290 .0(216.0) xX x 
¥2(D) 1,2,4,6 0.24 293 .0¢217 0) xX Xx 
¥3(M) 1,2,4.6 0.27 298 ,0( 267.0) x x 
R2 1 ? 300.0 x x 
P(B) 2 0.0 300.0 x x 
¥1(¥) 1,2,4,6 0.26 305 .0(273.0) x x 
LAL 1,2 0.55 310.0 x Xx 
OR3(D) 1,2,3,4 0.36 327 .0(315.0) x x 
OR2(C) 1,2,3,4 0.36 339 .0(330.0) -- -- 
ORL(A) 1,2,3,4 0.36 365.0(350.0) X x 
NJ (1) 2 0,0 410.0 x x 
RL 1 2 450.0 k x 
B 2 Q.1 550.0 XK x 
P(C) 2 0.0 600.0 x x 
LA3 1,2 8.8 1114.0 x x 
I 2 0.0 1200,0 x x 
LALL(K) 7 0.25 4500.0 Xx x 
RI(P) ? 0.33 8800.0 Xx Xx 


8Case nomenclature relates co that reported in the literature; numbers 
and/or letters that may be parenthetical following the geographical 
location keys given below, designate specific individuals. 


LA: Los Alamos A: Argonne RI: Rhode Island P: Picesburgh 
NJ: New Jersey R: Russia UT: U. of Tennessee OR; Oak Ridge 
B: Belgium T: Icaly C: Czechoslovakia Y: Yugoslavia 


Okey to References: 

Thoma and Wald, 1959. 
Hibner and Fry, 1980. 
Andrews et al., 1959, 
Braverstock and Ash, 1983. 
Klener ec al., 1986. 
Lushbaugh, 1969. 

Fanger and Lushbaugh, 1967, 


WOU wenre 


CM{dline body or mean bone marrow dose neutron (RBE = L). 
4nose estimates in parenthesis from Braverstock and Ash, 1983. 


“average dose to stomach and intestines. 
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B.2 ANALYSIS. 


The data in Table 31 were assumed to result from a cumulative 


distribution, 6(D; a,f), that represents the fraction of individuals 


with symptoms or signs after an acute dose. D. The two adjustable 


parameters, a and f, were determined by the maximum likelihood tech- 


nique [Cox, 1983], which uses the ¢:ta in binary form. We were then 


able to avoid the disadvantages of having to divide the data into 


arbitrary groups, as required for a standard regression analysis. 


For each of the 40 data points from Table 31, 


the fractional 


incidence, yj, for the ith point is unity if the symptom occurred (X), 


or zero if not (--). If the corresponding dose is 
likelihood function has a factor 6(D;) whenever yj 


of 1 - @(Dj) otherwise, i.e., 


Y; (1 - y,) 
L = I] a(D,) ~ . (1 - 6(D,)) 


where the product extends over the 40 points. The 
distribution @(D) then results from minimizing the 
this function, i.e., we find the parameters of the 


which minimize 


Dj, then the 


= 1, and a factor 


, (B.1) 


best fit to the 
negative log of 


distribution models 


-In L= - »: (Y; 1n@(D, ) + (1 - ¥) In(l - 8(D;))) . (B.2) 


Values for the two adjustable parameters, a and §, were ceter- 


mined Zor each of the five distributions considered, using a simplex 


algorithm [Press et al., 1986]. Table 32 gives the values for the a 


and § parameters. The log likelihood function was then used to develop 


and calculate the variance-covariance matrix for the optimized values 


of a and #. The resulting parameter standard deviations (0, and og) 


and covariances (cov(a, $)) are also included in Table 32, as are the 


functional forms for statistical models, @(D; a, 8). 


Table 32. Regression Models and Parameters 


Model a Ga B op cov(a,B) 
Normal -2.196 0.6798 0.01243 0.003220 -0.001946 
Set I Log-normal -11.944 4.0604 2.362 0.7654 -3.097 
Nausea Weibull -11.369 3.6275 2.158 0.6584 -2.380 
Logistic -4.055 1.4077 0.02391 0.007621 -0.009753 
Log-logit ~-22.473 8.6452 4.452 1.661 -14.32 


Normal -2.289 0.7089 0.01207 0.003206 -0.002053 
Set I Log-normal -12.024 4.1769 2.342 0.7788 -3.243 
Vomiting Weibull ~11.989 893.9672 2.243 0.7166 ~2.835 
Logistic -4.013 1.3610 0.02180 0.006751 -0.008375 
Log-logit -21.596 8.1033 4.214 1.534 -12.40 
Normal ~2.097 0.6575 0.01258 0.003326 -0.001946 
Set II Log-normal -11.363 3.8822 2.274 0.7407 -2.865 
Nausea Weibull -10.671 3.4010 2.048 0.6256 -2.120 
Logistic -3.988 1.4194 0.02487 0.008059 -0.01051 
Log-logit -21.352 8.2970 4.276 1.607 -13.30 
Normal ~2.236 0.7013 0.01248 0.00339 -0.002153 
Set I] _Log-normal -12.050 4.2532 2.369 0.8016 -3.400 
Vomiting Weibull -11.498 3.7530 2.175 0.6864 -2.568 
Logistic ~4.071 1.4277 0.02359 0.007621 -0.01006 
Log-logit -22.210 8.6266 4.374 1.651 -14,21 
4Models: 
Normal 8(D) = {a + BD) 
Log-normal 6(inD) = O(a + BlnD) 
Weibull 6(1nD) = 1 - exp {-e% + AlnD) 
Logistic 6(D) = 1/{1 + exp [-(a + fD)}}) 


Log-logistic @(1nD) = 1/{1 + exp [{-(a + BlnD)]}) 


B.3 RESULTS. 

Using the regression relationship, dose response calculations 
were performed based on the dual data sets (data sets I and II). Data 
set II differs from set I only for the dose values in parentheses in 
Table 31, reflecting the estimates of Baverstock and Ash [1983]. Dose 
response probability for nausea and vomiting are expressed in 
Tables 33 and 34. 

Table 33 summarizes the effective 10, 50, and 90 percentile doses 
based on data set I for the incidence of nausea and vomiting for each 
of the five statistical models. Also included are the 90 percent 
confidence limits for each value. Table 34 gives analogous informa- 
tion for data set II. Based on the x2 goodness-of-fit statistic, the 
normal and logistic distribution models seem to provide the least best 
fits of the data which is not surprising since those models do not 
predict zero incidence at zero dose. Actually, all the models fit 
both sets of data reasonably well, and which one that is chosen to 
represent the response becomes a matter of preference. 

Nausea and vomiting response curves are plotted in Figs. 92 and 
93 based on the analysis of data set I using the log-normal distribu- 
tion model (graphic results for data set II are not plotted due to the 
similarity with data set I). Doses in Fig. 92 are midline tissue (MLT) 
absorbed dose values, and in Fig. 93 they are free-in-air (FIA) ex- 
posure values. For convenience, dual plots are shown: the top row 
are plots of incidence with linear probability (vertical) scales, 
while the bottom row has nonlinear (probability) scales that produce 
straight line plots for incidence; in both cases, the abscissas are 
log dose. The 90 percent confidence bounds (dashed lines) are also 
indicated in the plots. 

The question of neutron RBE still remains an open one, and it is 
unlikely that it will be resolved empirically for humans short of 
possibly future clinical experience utilizing neutron radiation 
therapy. With regard to laboratory work with monkeys, Young [1986] 
has suggested that neutrons are more effective than gamma radiation in 
producing emesis. On the other hand, based on reactor irradiation 


studies with dogs, Cordts et al. [1985] found neutrons less effective 
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(RBE = 0.48) in producing emesis. Before a more definitive applica- 
tion of RBE can be specified for emesis in humans, further animal 
experimentation will be required together with a more thorough 
deiineation of the radiation-induced mechanisms. 

As an exercise to illustrate the effect of neutron RBE, we per- 
formed some additional likelihood calculations utilizing a lognormal 
model to analyze the accident data (data set 1). Based on the n/y 
ratios given in column 3 of Table 31, we assumed various values for 
the neutron RBE ranging from 0.5 to 2.0 and reanalyzed the vomiting 
response for those conditions. The calculated neutron RBE effect 
obtained is shown in Fig. 94 for the EDsg emesis endpoint both in 
terms of MLT and FIA dose; it amounts to about a 33 percent increase 


in the EDsg for a factor of four decrease in RBE. 
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Figure 94, EDsg for emesis versus neutron RBE (basis: data set 13 
lognormal] model). 
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APPENDIX ¢ 
EQUATIONS OF THE GUT INJURY MODEL 


This appendix presents the equations for the gut injury model 
(GIM). The physiological and anatomical basis for the model and the 
structure of the model are discussed in Sec. 5 of this report. 

The foundation for the radiation response of the GIM is the 
lethal potentially lethal (LPL) model. Curtis [1986] has published a 
full discussion of the LPL model and its application to a wide range 
of radiobiological response data. That model includes two types of 
radiation-induced chromosome lesions, lethal lesions, and potentially 
lethal lesions. The lethal lesions, having an average number of no 
per cell, are irreparable and prevent mitosis. A cell may have more 
than one lethal lesion, but any one is sufficient to cause reproduc- 
tive death of the cell. The potentially lethal (PL) lesions, having 
an average number of ng per cell, are repairable by a process based on 
linear kinetics, All PL lesions can be repaired in time; however, a 
cell that enters mitosis with an unrepaired PL lesion will fail to 
divide and be reproductively dead. The LPL model also includes a 
misrepair term that generates lethal (irreparable) lesions from PL 
lesions at a rate that is quadratic in the concentration of PL 
lesions. This nonlinear term produces a shoulder on the cell survival 
curve and provides a link to the linear-quadratic (LQ), or alpha/beta, 
model frequently used to analyze cell survival data. 


The equations for the LPL model are: 


R-¢«. n -e n : (1) 


Se 
tT 


R +e n i (2) 


where R is the dose rate. Both types of lesions are produced at a 
rate proportional to R. Figure 95 provides a diagrammatic definition 


of the rate coefficients and variables in these equations and shows 
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(No damage) Ng = Mean number repairable iesions/cell 


So = mean number irreparabie lesions/cell 


x 


TBA frepair) LPL model LO modet 
\ 
(Repairable RAB ‘ Nac a 
damage) 
2 
n 
A 
NaC : 8 
(rreparable 2f€ an - Ene) 
damage) BA "BC 
és F BA u 
BC (misrepair) 
4 ¢ 
c ese 


Cell survival fraction: S(D. t} = exp (-Mg NG) 


Source: Thames (1985), based on Tobias et.al. 41980), and Curls (1982). 


Figure 95, Cell radiation resronse based on the lethal, 
potentially lethal (LPL) cell lesions model. 


the correspondence with the parameters customarily used in the LQ 
model. 

In the LPL model, the distribution of each type of lesion among a 
group of cells is assumed to follow a Poisson distribution since the 
lesions are produced at random by radiation exposure. This assumption 
leads to the exponential expression for cell survival shown in 
Fig. 95. The expression is simply the probability that a given cell 
has no lesions of either kind at the time it enters mitosis. 

We had to address a major problem crucial to applying the LPL 
model to a proliferating tissue. Previously, the LPL model had been 
applied only to data where the exposure and repair times were short 
compared to the cell cycle. The difficulty for longer times arises 
from the fact that proliferation disrupts the Poisson distribution of 
lesions. Successful cell division occurs in cells with no lesions, 
lowering the average number of lesions per cell but not in a random 
manner. Likewise, the death of a cell entering into mitosis due to 
the presence of one or more lesions removes any or all lesions 
present. We therefore had to develop equations that would track both 


lesion number and cell number over long periods of time given any 
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exposure history. In addition, we needed to replace the simple cell 
survival formula shown in Fig. 95 with a formulation based on these 
same principles. 

Probably the most straightforward analytic approach to this 
problem would be a Monte Carlo procedure for tracking a large number 
of cells. The distribution of lesions would then evolve into a non- 
Poisson form as it must do so in reality. To avoid the extensive 
calculation requirements of a Monte Carlo model, we derived an ap- 
proximate solution based on continuous time differential equations. We 
refer to our formulation as the proliferation and intracellular repair 
(PAIR) model. 

The key to this approximate solution was to divide the 
proliferating cells into three classes, labeled A, B, and C, ina 
manner similar to Curtis's classes of lesions. However, it is impor- 
tant to point out that in the PAIR model, the classes refer to cells 
and not to lesions. The number of cells in each class and the defini- 


tions of the classes are as follows: 


Na = number of cells with no lesions of either kind (uninjured 
cells), 

Np = number of cells with one or more PL lesions, but no lethal 
lesion (injured cells), and 

Nc = number of cells with at least one lethal lesion (mitotically 
dead cells). 


Proliferation increases the number of cells in class A in the normal 
manner but does not occur in class C since those cells are mitotically 
dead. Proliferation (attempted division) of class B results in cell 
death and moves them to class C since we assume, following Curtis, 
that any cell entering mitosis with a PL lesion suffers mitotic death. 
It is not necessary to track the number of lesions in class C 
cells since they are mitotically dead. However, we track the number 
of C cells because they are still present in the tissue. We assume 
that their presence affects homeostasis and that they differentiate to 


form mature, functional cells. These assumptions seem to fit the 
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dynamics of the intestinal epithelium but may not be appropriate for 
every type of tissue. 
The only lesion number that must be tracked in the PAIR model is 


that of the class B cells, where: 
ng = mean number of PL lesions in class B cells. 


Note that by the definition of class B cells, ng is always greater 
than or equal to one. Therefore, ng is not the average number of PL 
lesions as defined in the LPL model since that number goes to zero in 
the absence of radiation exposure. For convenience, we define a 


different number that does correspond to the LPL model: 


Npl = mean number of PL lesions in a hypothetical pool of cells 
with a Poisson distribution of PL lesions and for which the 
number of cells with one or more PL lesions and no lethal 


lesions is Np. 


In other words, we approximate the actual distribution of lesions in 
the B class of cells by assuming that it is a truncated Poisson dis- 
tribution from this hypothetical pool of cells. The accuracy of this 
assumption remains to be investigated. The advantage of the assump- 
tion is that proliferation of class A cells can occur without upset- 
ting the statistical distribution of lesions in class B cells. Fur- 
thermore, in the limit of no proliferation, the mathematical results 
match those of the LPL model. 


The value of np in terms of npy is given by: 


n= ——— (3) 


The rate equations for the average number of PL lesions per cell 


in the injured compartment is: 
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* Na 
Np, = MapRE - ("3 = 1) Ne MapRE 


Exposure Dilution 
(4) 
~ €Ba™PL ~ €gcpr"s(2 i ("5 = Mp) |E ; 
Repair Misrepair 
where R is the radiation exposure rate (Gy/h), 
PL 
eG ny 1 - [n, - 1 y 
B PL 
(3) 
and lim E = lim =2 . 
ny (ny 1/2) 
Np 70 nl 
The rate equation for the number of uninjured celis Ng is: 
Ny =  ANg = (M3 * Mac} Rg 
Proliferation Exposure Losses 
(6) 
+ €,, ("5 - Np) Np z ON, 


Repair of Injured To Pipeline 
Cells 


The coefficient 6 which determines the rate at which clonogens feed 


the "pipeline" to the transit compartment is defined below, as is the 
pip iP 


clonogen division rate i. 


The rate equation for the number of injured cells Np (potentially 
lethal injury only) is: 
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Ma ANG! > ag RNs MacRNp 
Mitotic Exposure of Exposure Causing 
Death Uninjured Cells Lethal Injury 
(7) 
- pa("a - "pr)%a ~ Sec™pr™aNn > ONG 
Repair Misrepair Pipeline 


The rate equation for the number of killed cells No (lethally 


injured or mitotically dead, but still functional) is: 


we 


c= MacR(Nq + N3] + €, 0M) TRN, + AN, - ON. (8) 


Exposure Misrepair Mitotic Death of Pipeline 
Injured Cells 


The instantaneous clonogen division rate \ is affected by both 


homeostasis and mitotic delay as given by: 


E, 


(9) 


where, 
H = homeostasis factor(dimensionless) 
M = mitotic delay factor (dimensionless) 
y = normal division rate to balance attrition (holy, 


Normal conditions give H = M = 1 and A = y. The maximum value of 
X is Am, determined by H. The minimum value of \ is 0, 
determined by M. 

Radiation exposure slows down cell division through a mitotic 
delay process. For the GIM, we have modeled this process based on a 
saturable enzyme repair concept where a hypothetical damage level Q is 
reduced by the action of a finite pool of repair enzymes. Assuming 
cell cycling proceeds normally when the repair enzymes are not 


activated, and slows down wien they are activated, we employed a 
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Michaelis-Menten form to describe activated enzyme saturation. This 


results in our rate equation to express cellular damage piven by: 


° a Q 
Qi= R ° 5) AQ (10) 

Damage Damage 

Production Repair 


Q = cellular damage (Gy) 

Am = freely growing (maximum) cell division rate (h71) 
6 = acute dose mitotic delay constant (Gy) 

A = threshold for saturation of repair = 0.1 Gy 


R = dose rate (Gy/h) 
The mitotic delay factor M in Eq. (9) above is then given by: 
M = [1-tanh(Q/A)] (11) 


When the damage level Q is low, M approaches one and cell division 
proceeds normally. When the damage level increases beyond the charac- 
teristic value A, M rapidly approaches zero according to the hyper- 
bolic tangent function, effectively halting the progression of cell 
cycling. 

An H-tissue such as the intestinal epithelium must have com- 
munication mechanisms between compartments in order to maintain 
homeostatic equilibrium. Figure 81 in Sec. 5 illustrates the control 
mechanisms that govern the GIM model. Many links are possible, and 
the nature and number of links have not yet been established 
experimentally. We have attempted to choose the simplest and minimum 
number of (feedback) control links that will both maintain equilibrium 
in the three GIM compartments and provide a correct qualitative 
response to radiation insult. 

The primary homeostasis loop is the control of clonogen 


proliferation through the division rate \. When the tissue is in 


equilibrium with Ng clonogenic cells, the normal division rate is y 
expressed in units of divisions per hour per cell. Compensatory 
proliferation is accomplished through the homeostatic control 
factor H. In equilibrium, both factors M and H are equal to l, and 
the production rate of new clonogens is Nj, just balancing the attri- 
tion rate from the villi to the intestinal lumen. 

The mathematical form of the homeostasis factor H is presented in 
Fig. 96. The primary purpose of H in the GIM is to increase the cell 
division rate to its maximum value Aq when the population of the 
transit compartment drops to zero. At first, we used a linear depend- 
ence on Nr, but found that the delayed recognition of damage was not 
properly represented. This phenomenon shows up in both fractionated 
[Withers and Elkind, 1969] and constant dose rate (Withers, 1972] 
studies. In spite of cell depletion, it is found that rapid 
proliferation is delayed for at least two days. We introduce this 
behavior by using a small value for the exponent a = 0.21, to moderate 
the increase in cell division rate as Ny decreases. A final 
factor £(N/No), shown in Fig. 96, is used to limit the growth of the 
clonogenic compartment to twice its normal value. As explained later, 


Ny is not allowed to go above its equilibrium value Nz g, so that the 
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Figure 96. Clonocen homeostasis factor, H, 
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factor H will not drop below 1. This restriction ensures a supply of 
cells to balance the normal attrition rate unless overridden by 
mitotic delay. 

The GIM medel uses the transit population Ny rather than the 
villi population Ny to control compensatory proliferation mainly on 
the grounds of physical proximity. We do not have experimental data 
to support this choice. In fact, it might be useful to include Ny in 
the equation for H as an attempt to improve the description of delayed 
recognition of damage. 

The remaining control mechanisms shown in Fig. 81 of Sec. 5 
involve the fluxes of cells in and out of the three compartments. 
Figure 97 shows the mathematical expressions for the fluxes. The 
attrition rate from the villus has the simplest mathematical form. It 
is held fixed at the normal level yNg when the villus population is 
between 0 and its equilibrium value Ny,. When the population is zero, 
the villus output is set equal to its input if the input is less than 
yNo and is equal to yNg otherwise. The assumption is that the villus 
has atrophied completely and will not grow until the supply of cells 


from the transit compartment exceeds the normal attrition rate. When 
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Figure 97. Cell fluxes between compartments in the GIM model; note 
that ¢ is subscripted, but 8 is not in the equations. 
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the population reaches its full value Nyg, any input level that ex- 


ceeds yNo is passed on through in order to limit the villi to their 
normal size. 

The flux of cells from the transit compartment to the villus is 
increased linearly from its normal value yNg to twice that value as 
the villus level drops from normal to zero. Like the villus, the 
transit compartment is assumed to shrink to zero size when it is empty 
so that output is equal to the input if the input is less than the 
output demanded by the villus based on the function f"(Ny/Nyg). Also, 
when the transit compartment reaches full size, it passes excessive 
inputs on through in order to limit its size to Ny. 

The output from the clonogen compartment is more complex. Like 
the others, it is designed to provide a flux equal to the attrition 
rate yN, at equilibrium. The flux @N shown in Fig. 81 of Sec. 5 is 
defined in Fig. 97. The mitotic delay factor M appears in the output 
flux @6N as well as in the division rate 4. The assumption is that if 
mitotic delay turns off cell division, it will also turn off cell 
differentiation that moves cells from the clonogen compartment to the 
transit compartment. Otherwise, mitotic delay in the model causes an 
anomalous depletion of the clonogen compartment that is not consistent 


with experimental data. 
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